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ABSTRACT 

Context. The study of high redshift Tully-Fisher relations (TFRs) is limited by the use of long slit spectrographs, rest frame B band 
and star formation selected galaxies. 

Aims. We try to circumvent these issues by using integral field spectroscopy (SINFONI), by studying the rest frame K band and stellar 
mass TFR, and by selecting targets without a bias to strongly star forming galaxies. In this paper, we demonstrate our methods on our 
best case. This galaxy, F257, at z = 2.03, was selecte from a sample of candidate high redshift large disk galaxies in the Hubble Deep 
Field South that were selected with photometric and morphological criteria. 

Methods. We used SINFONI at the VLT to obtain an integral field spectrum of the Ha line and hence a velocity field and rotation 
curve. We also use UBVIJHK+IRAC band photometry to determine a stellar photometric mass. 

Results. We find that F257 is indistinguishable from local late type galaxies in many respects: it has a regular velocity field, increasing 
velocity disperion towards its center, its rotation curve flattens at 1-2 disk scale lengths, it has the same specific angular momentum 
as local disks, its properties are consistent with the local K band TFR. Although mainly rotationally supported, its gas component is 
dynamically heated with respect to local galaxies (V/cr. ~ 4) and it is offset from the local stellar mass TFR at the 2tr level. But, this 
offset depends on the SED modeling parameters. In particular, for a 2-component star formation history (SFH), F257 is in agreement 
with the local stellar mass TFR. F257 is then a nearly (~ 75%) maximum disk. The dynamical properties of F257 are more like those 
of local galaxies than those of any other galaxy at similar redshift observed to date. However, the gas-to-stellar mass ratio is unusally 
large: 2.5. 

Key words. Galaxies: high-redshift - Galaxies: kinematics and dynamics - Galaxies: evolution 



m ' 1 . Introduction More than 30 years after the first publication of what is now 

K> ' . . , , . . , , . , , known as the Tully-Fisher relation (TFR) (Tully & Fisher 1977 ), 

, Starting with density fluctuations in the early universe, dark mat- ±e Q1 • ifl of ^ TFR k &m subject of debate Some authors ex _ 

g . ter halos grow through gravity. Galaxy formation is described by Mn ^ ppR af , a con nce of se lf. r eg U lated star formation 

semi-analytical models or hydrodynarmcal simulations. These fa difjks ^ masses ( snk rj997}. i n other mod els, 

may include physical processes such as gas cooling star for- ^ ppR h a direct co ence of the cosmo i ogica i equivalence 

mation, feedback from ste ar winds, supernovae and/or AGN between masfj and drcular velod resuhi from the finke 

chemical evolution and stellar populations (Okamoto [2007] and of ^ uniyerse which imposes a maximum radius from with mat . 

references therein). ter can accrete to form a disk (e.g. Mo et al. 1998, Steinmetz & 

The disk formation epoch in earlier generations of models Navamj m9 ) Steinmetz & Navarro ( 1999 j show mat akh h 

is relatively late and mostly at z S 1. These early models a so ^ d rf ^ ppR c£m bg lained b hierarc hi- 

predicted that these disks are much smaller than observed in lo- ^ modeK ^ normalization and evolution of the TFR depend 

cal spiral galaxies and that their specific angular momenta are j on the star formation preS criptions and the cosmologi- 

only 20% of that of their halo's (Navarro, Frenk & WluteM- cal parame ters. Recently, Courteau et al. showed that two 

More recent models produce more realistic size disks (Mo etal j knQ ^ different theoretical pred i ctions for two different 

1998 Portinari & Sommer-Larsen|20El Governato et al.|2SD slopes f the TFR are in fact related, 
and place the formation epoch somewhat earlier: 0.4 < z < 2.0 

(Buchalter, Jimenez & Kamionkowski 2001 ). Reproducing the Dlsk § alax y scalln g relations, like the TFR, put strong con- 
specific angular momentum of local disk galaxies remains a straints on § alax y formation scenarios, and in particular on the 
challenge (e.g. Portinari & Sommer-Larsen|2007). connection of the dark and baryomc component (van den Bosch 
2000). The zero-point of the TFR constrains the SFH and the 



Send offprint requests to: cosmoiogicai parameters that control the universal baryon frac- 

* Observations carried out using the Very Large Telescope at the ESO tion and the time of assem bly of galaxies of different mass 
Paranal Observatory, Program ID 077.B-0079. (Steinmetz & Navarro 119991 ). Many authors struggle to repro- 



2 



L. van Starkenburg et al.: Dynamical properties of a large young disk galaxy at z=2.03 



duce the zero-point of the TFR (e.g. Portinari & Sommer-Larsen 
120071 Steinmetz & Navarro [T999j l. The slope of the TFR and 
the steepening of the slope towards longer wavelengths can only 
be explained by supernova feedback (van den Bosch 200^1 al- 
though others claim it is a natural result of the hierarchical mod- 
els (Steinmetz & Navarro [l999l ). The small scatter in the TFR is 
explained by scatter in halo spin parameters and formation red- 
shifts (e.g. van den Bosch 2000). Even the most recent galaxy 
formation and evolution models are not able to simultaneously 
reproduce for example the slope, zero-point and scatter of the 
TFR and the zero-point and shape of the luminosity function 
(e.g. Courteau et al. 2007 and references therein). 

The TFR defines a tight correlation between the luminosity 
and the fiat rotation curve velocity of spiral galaxies. If luminos- 
ity is converted into stellar mass (which requires detailed mod- 
elling of the stellar spectral energy distribution (SED)), the tight 
correlation remains: the stellar mass TFR. Semi-analytic models 
predict evolution in the slope, zero-point and scatter properties 
of the TFR and the stellar mass TFR (e.g. Steinmetz & Navarro 
1999 Buchalter, Jimenez & Kamionkowski 2001 Portinari & 
Sommer-Larsen 2007). A key element is the relative evolution 
of the dark and stellar mass components, which is a direct test of 
the stellar mass buildup within a dark halo of a given mass. 

Critical tests of these theories can be made by observing the 
evolution of disks backward in time. Rotation curves (RCs) of 
field galaxies out to z = 1 measured with the Keck telescope and 
Very Large Telescope (VLT) are found to be similar to those of 
nearby galaxies in both shape and amplitude implying that some 
massive disks were already in place at z ~ 1 (Vogt et al. 1996, 
[19971 Barden et al. 120031 Swinbank et al. l2006l l. Observations of 
z ~ 2 galaxies have revealed ordered rotation, but only the rising 
part of the rotation curve (e.g. the lensed sources of Lemoine- 
Busserolle et al. 2003 ). Van Starkenburg et al. (2006) found or- 
dered rotation in a z = 2.46 galaxy, with an asymmetric line 
profile indicating a central starburst and weaker starformation in 
the outer parts. Ordered rotation has been found out to z ~ 3.2 
by Moorwood et al. d2000l [2003b and by Nesvadba et al. ( [20061 ), 
the latter again on a lensed source. 

Also, several attempts to measure size evolution and the TFR 
at redshift 0.5-1.5 have been made over the past few years with 
contradictory claims about disk evolution. For instance, Mo et 
al. dl 9981 l use d prev iously acqui red da ta (Vogt et al l 19961 [19971 
Schade et al. 119961 Lilly et al. (1998) to conclude that z ~ 1 
galaxies are smaller than local galaxies at a fixed rest-frame B- 
magnitude or rotation velocity, consistent with hierarchical mod- 
els. A brightening with redshift of disk galaxies at fixed rota- 
tion speeds was also detected by Barden et al. (12003). However, 
Simard et al. (1999 ) and Ziegler et al. (120021) have shown that 
massive galaxies at z < 1 are consistent with having no evolu- 
tion in surface brightness, once selection effects are taken into 
account. 

The largest source of uncertainty in those studies, all done 
with long slit spectra, are the velocity measurements. In recent 
years, it has become possible to do integral field spectroscopy 
at near infrared wavelengths, for example with SINFONI at 
the VLT (Eisenhauer et al. |2"U03l Bonnet et al. 120041 ). Integral 
field spectroscopy has significant advantages over long slit spec- 
troscopy. First, long slit measurements may suffer from slit mis- 
alignment. More fundamentally, the photometric and kinematic 
axes may be different (due to star forming regions or intrinsic 
shape of galaxy). Kinematic disturbances are harder or impossi- 
ble to identify in slit spectra. Integral field spectroscopy does 
not have all these slit related problems. Finally, integral field 
spectroscopy makes for a fairer comparison with low-z Tully- 



Fisher data, which are typically derived from full imaging of 
the HI velocity field with radio aperture synthesis. Integral field 
spectroscopy thus presents a major improvement over slit spec- 
troscopy, provided significant resolution over the rotating disk is 
achieved. Recently, Forster Schreiber et al. (2006) and Bouche 
et al. (2007 ) used SINFONI to obtain two-dimensional velocity 
fields of a number of galaxies at z ~ 2 with different selection 
criteria, clearly demonstrating the power of an efficient integral- 
field spectrograph such as SINFONI for this type of work. 

A second major problem that plagued previous high-z TFR 
studies is the use of rest-frame B-band, which is very sensitive 
to small amounts of recent star formation and extinction by dust. 
There are now a number of deep surveys available with compre- 
hensive photometric coverage ideal for studying the high redshift 
TFR in rest frame near-IR luminosities. We used the HDF-S and 
MS 1054 fields which includes recent Spitzer IRAC data (Labbe 
et al. 120051 ). This allows us to use rest-frame near-IR luminosi- 
ties, which define much tighter TFRs than rest-frame optical lu- 
minosities (e.g. Verheiien l2001l) . Rest-frame K band luminosi- 
ties as derived from Spitzer data are much less sensitive to the 
effects of dust and recent starformation than rest frame optical or 
UV luminosities and will provide us with more accurate stellar 
masses. 

A third limitation of previous studies is that many of these 
target starformation-selected galaxies (e.g., selected on rest- 
frame UV or bright emission lines). This compromises attempts 
to use these galaxies to trace the evolution of the TFR (see 
van Starkenburg et al. 2006 for a detailed discussion), since lo- 
cal TFR samples are composed of relatively quiescent galaxies. 
Although a certain amount of star formation is necessary to ob- 
serve the velocity fields of high redshift galaxies, there are sev- 
eral ways to partly circumvent this problem. 

First, by selecting from ultra-deep near-IR imaging surveys, 
such as the FIRES data of the HDF-S and MS 1054 fieldsQUabbe 
et al. 120031 Franx et al. 120031 Forster-Schreiber et al. 20061 ), we 
can select galaxies in the rest frame optical instead of the rest 
frame UV. The FIRES data is sufficiently deep to select galaxies 
in the rest-frame optical out to z ~ 4 and the depth on the Hubble 
Deep Field South (HDF-S) is sufficient to detect non-evolving 
(0.6L, at z — 0) spirals out to z = 3. 

A second improvement is to preselect massive galaxies. 
Contrary to low mass galaxies, their log(L) (and hence offset 
from the TFR) is not or less affected by the luminosity of a star- 
burst (van Starkenburg et al. 2006). We observed several targets 
in this category in two redshift ranges, which will be discussed 
in two separate papers (van Starkenburg et al. 2008a, 2008b). 

Finally, it is now possible to select high redshift galaxies that 
morphologically resemble local disks. A search for such disks 
in the FIRE survey of the 2.'5 x 2'.5 HDF-S field resulted in a 
surprisingly large number (6) of candidate large disks (Labbe et 
al. 120031) . These galaxies, at redshifts 1.4 to 3.0 (most of them 
spectroscopically confirmed), have exponential profiles with ef- 
fective radii ranging from 0765 to 079 (or 5.0 - l.Sh^ kpc) and 
extending over 2-3 effective radii and appear remarkably regular 
in the near infrared (rest frame optical). Four are well fit by just 
exponential disks and two require an additional r 1 ^ 4 bulge com- 
ponent, contributing about 40% of the light. At shorter wave- 
lengths, they have more extended and irregular morphologies. 
They have strong rest-frame optical breaks in their centres and 
evidence for intense star formation in the outer parts, reminis- 
cent of local spirals, with blue star-forming disks and red bulges. 



1 The FIRES data are publicly available, see 
http://www.strw.leidenuniv.nl/~fires 
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Their sizes are comparable to the Milky Way and those of L* 
disks in the local universe. Specifically, these sizes are much 
larger than those of typical LBGs (1 - 2h^ Q l kpc; Giavalisco et 
al.[l996; Lowenthal et al. 1 19971 )) and their sizes are in the upper 
end of the range found by Trujillo et al. d20061 l for all galax- 
ies at similar z in the two FIRES fields. They are bright in K 
(K= 19.7 - 21.1) and have large stellar masses (M > 10 11 M G ). 
They constitute half of the most rest-frame luminous galax- 
ies (Ly > 6 x 10 10 'h^L Q ) and their number density exceeds 
model predictions (e.g., Mo et al., [1998 ) by a factor of two. 
Obviously, the first and foremost question that needs to be re- 
solved is whether these objects are truly rotationally supported 
disks, and this can only be tested by obtaining kinematic data. 

In this paper, we present SINFONI observations of one of 
the large disks of the Labbe et al. sample (id 257 in the FIRES 
HDF-S catalog). These observations were part of a larger project 
of SINFONI observations of mass and/or disk selected high red- 
shift galaxies, which will be presented in two separate papers 
(van Starkenburg et al. 2008a 2008b ). Here, we will discuss our 
methods as used on our best case. 

We observed F257 with SINFONI during the night of August 
15, 2006. We selected F257 as our primary target for this ob- 
serving run because it is the second brightest K band source 
(K = 20.25) of the large disks in the HDF-S (Labbe et al. 120031 ) 
that has a known spectroscopic redshift (z = 2.03, Vanzella et 
al. 120021 ) with Ha at an observable wavelength with SINFONI. 
(We had observed the brightest K band target fullfilling this cri- 
teria during an earlier observing run, see van Starkenburg et al. 
2008b). Its effective radius is C'74 in K band, the total radial ex- 
tent above the K-band detection threshold is two effective radii 
or 3-3.5 scale lengths (R d = 0744). The in clinati on (~ 54°, cal- 
culated from the ellipticity of Labbe et al. 2003 ) is suited for a 
velocity field analysis. F257 shows a clumpy distribution in the 
HST F814W filter (see Fig.QJ) which corresponds to rest frame 
UV. The rest frame optical (observed K band) image shows a 
smooth exponential profile (see Fig. QJ). This difference can- 
not be attributed to the different PSF of the F814W and the K 
band image: the F8 14W image smoothed to the K band PSF still 
shows the clumps. Labbe et al. ([2003 ) fitted an exponential pro- 
file without the need for a bulge component. NICMOS H band 
shows an intermediate picture (see Fig.QJi). The Sersic index as 
determined by Trujillo et al. ( [20061) o n the FIRES (ISAAC) H 
band image is 0.76. Wuyts et al. (2007 ) made an extensive study 
of the SEDs of high redshift galaxies in HDF-S and found for 
F257M, = 3.5+°^xl0 10 M G , a mass weighted age of - 160 Myr 
and Ay — 1.0 using a Salpeter initial mass function (IMF), solar 
metallicity, a t 300 model for the star formation history (SFH) 
and Bruzual & Chariot (2003 ) models. A number of z ~ 2 galax- 
ies have been observed with SINFONI and other IFUs (Genzel 
et al. 120061 For ster-Schreiber et al. [20061 Bouche et al. [20071 
Law et al. 120071 ). F257 is not different in terms of colors, SFR or 
age, and satisfies the conventional U-dropout criterium. It is the 
only galaxy selected as large disk galaxy, but the other galax- 
ies were selected from catalogs with insufficient resolution to 
identify morphological large disk galaxies at optical rest frame 
wavelengths. 

Throughout this paper, we use Ho = 70 km s~'Mpc~' , fl^ = 
0.3, Qa = 0.7 and Vega magnitudes. 




Fig.l. From top to bottom: F814W image, NICMOS H band 
image and Ks image from FIRES data. Image size is always 
3"times 4"and the images have been rotated so that east is to 
the top and north is to the right. The (K band) center from the 
FIRES data is marked by the white crosses. 



2. Observing strategy 

SINFONI (Eisenhauer et al. 120031 Bonnet et al. |2004l ) is the near 
infrared integral field spectrograph of the VLT. SINFONI slices 
a square image into 32 slitlets that are combined into a single 



pseudo slit and then dispersed. It operates in J, H, K and H+K 
band and offers three image sizes. We used the 8 x 8" field of 
view (FOV) in K band. The spectral resolution in this mode is 
R ~ 4500 or FWHM ~ 75 km s" 1 . 
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We used an ABA' B' -like observing cycle, comparable to the 
ABBA-cycles used in near-IR long slit spectroscopy. Placing our 
targets at different positions in the FOV instead of taking sky 
frames ensures optimal use of observing time. Because there 
is a pair of galaxies close to F257, the last two offset posi- 
tions are only slightly offset from the first two positions. We 
rotated SINFONI's FOV so that the other galaxies do not in- 
terfere with the observations of F257 (see Appendix [A] for de- 
tails on SINFONI's FOV). Individual exposures were 900 s, to- 
tal on source integration time 6 hours. We used 900 s exposures 
to keep the extra overheads (observations of PSF stars, see de- 
tails in Appendix [Ab within reasonable limits. As our target is 
too faint to acquire directly, we performed a blind offset from 
a bright nearby star to our target. The seeing varied during the 
night between 0735 and 0786 (measured on the PSF stars). The 
seeing on the weighted combined cube is 0755. We also observed 
four hot standard stars for flux calibration during the entire night 
(which also included observation of another target). 

3. Data reduction 

We used a combination of IRAF tools and the SINFONI pipeline 
(version 1.3.0) (Modigliani et al. 2007) for data reduction. We 
reduced the standard calibration files provided by ESO for 
SINFONI (dark, distortion, flat, linearity, and arc frames) us- 
ing the SINFONI pipeline to get bad pixel maps, darks, flats, 
wavelength calibration map, the positions of the slitlets on the 
detector and distortion correction. 

We removed bad pixels from the science images and used 
LA_COSMIC (van Dokkum 2001) to remove cosmic rays. Some 
cosmic rays close to bright OH lines remained, we removed 
those close to the Ha line manually by replacing those pixels 
with the interpolated value of the surrounding pixels. We noticed 
an odd-even effect in some slitlets, with strength varying with 
time and in different quarters of the detector, sometimes absent. 
Therefore, we averaged all columns with the mean of their right 
and left neighbour to remove the odd-even effect. This causes 
a small (5% increase in the FWHM of the final PSF) smear in 
the x-direction after image reconstruction (see Appendix [Al for 
definition of x and y direction on SINFONI's FOV). 

We used the SINFONI pipeline to correct our data for distor- 
tion, flat fielding, reconstructing the data cube, and combining 
our science data in pairs ("first subtraction") and did exactly the 
same for the PSF stars. We did an illumination correction on the 
science data using an early version of the recipe made by Juha 
Reunanen. We checked the combined cubes for faint star-like 
persistence effects (see van Starkenburg et al. 2008a) but there 
were none. 

We measured the position and FWHM of the PSF stars at 
the wavelength range (averaged over 40 pixels) of the Ha and 
[Nn]/i6583 line. We used these positions to combine the data in 
pairs ("second subtraction") for reasons described in Appendix 
lAl The last two science frames were observed at the end of the 
night. There was no time left to observe PSF stars for these 
frames. We used the positions of the previous observation block 
(OB) for this pair (which will get very low weight in the final 
combination anyway because the seeing was poor). 

Then, we flux calibrated the data using the PSF and stan- 
dard star data. The flux of the PSF stars correlates well with 
1/airmass, and we used this relation to do a relative flux calibra- 
tion for all pairs of science frames. We also observed four bright 
standard stars during the night, but those did not show this tight 
correlation between airmass and flux (after correcting for differ- 
ent integration times and luminosity). Two standard stars gave 
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Fig. 2. The integrated ID spectrum of F257. We detect 
Ha and [Nn]^6583. The positions of the [Nn]/16548 and 
[S n]/16718/6733 are also marked in grey. The logarithm of the 
sky spectrum with OH lines is shown in the bottom panel. 

consistent results (with the PSF star relation), the third differed 
in the amount of atmospheric absorption, the fourth was more 
noisy. Fortunately, Ha is not at a relevant wavelength for atmo- 
spheric absorption. We choose to use the average of the two con- 
sistent standard stars for absolute flux calibration and ignored the 
other two standard stars. 

We used the FWHMs of the PSF stars to weigh the data when 
combining all data. We found the best result in terms of resolu- 
tion and signal-to-noise (S IN) when weighting with the square 
of the FWHM in the y direction (which is unaffected by the com- 
plex FOV reconstruction). We measured the seeing in the com- 
bined image from the combined PSF star image (using the same 
weights, at wavelength of Ha and [Nn]/16583) and found 0755. 
The FWHM in the x direction is 20% larger. 

Further details on the reconstruction of the FOV using the 
PSF stars are provided in Appendix lAl 

4. Results 

The final integrated ID spectrum of F257 is shown in Fig. 
[2] and shows Ha at redshift 2.028 together with [Nn]/16583 
([Nn]/16548 falls on top of a bright OH line). We do not detect 
the [S n] lines. The continuum is very faint but clearly detected 
in the median collapsed cube (excluding the wavelength range 
of the emission lines). The rest frame equivalent width (EW) of 
Ha is ~ 160 A where the continuum flux was calculated from a 
linear interpolation of the observed H and K magnitudes from 
the FIRES data. 

The [Nn]/Ha ratio is ~ 0.22 (averaged over the whole 
galaxy) which is consistent with ionization by star formation. 
The metallicity calculated from this ratio is 12 + log(0/H) = 8.6 
which corresponds to slightly subsolar metallicity (Denicolo 
et al. 120021 solar metallicity is 12 + log(0/H) = 8.91). The 
[Nn]/Ha ratio is slightly larger in the centre of the galaxy 
([Nn]/#a ~ 0.4) which gives metallicity 12 + \og{0/H) = 8.8. 
The [N n]/Ha ratio is also consistent with ionization by star for- 
mation in the centre of the galaxy. Ionization by an AGN would 
would give larger a larger [N n]/Ha ratio. 

Using the Kennicutt (1998) relation for a Salpeter IMF, 
the Ha flux corresponds to a star formation rate (SFR) of 
116M yr _I uncorrected for extinction. Using Ay = 1.0 from 
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Fig. 3. Synthetic long slit spectrum. Scales are indicated in the 
figure. Note the presence of faint [N n] emission. 
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Fig. 4. Rotation curve as measured from the synthetic long slit 
spectrum by fitting a gaussian to each spatial position in Fig. [3] 
The formal fit errors are comparable to the size of the symbols, 
except for an increase in the outer parts. The points are correlated 
due to the 0755 seeing. 



Wuyts et al. d2007b . the SFR becomes 232 M yr" 1 . This value is 
very close to the value found when dividing the best fit mass and 
age from the SED fitting: 219 M yr _1 . 

The FWHM of the Ha line is 244 km s after correction 
for instrumental resolution. Again, there is no evidence for the 
presence of an AGN in the form of broad lines. Stepping through 
the wavelength slices of the cube, one notices immediately the 
wavelength shift of He. A synthetic long slit spectrum (along the 
major axis, position angle 32° counterclockwise from the north) 
made from the cube integrated over the minor axis of the galaxy 
with slit width 1'.'25 shows the tilted emission line that flattens at 
at least one side, see Fig. [3] The RC measured from the synthetic 
long slit spectrum is shown in Fig. [4] We will discuss the full 
velocity field analysis in the next section. 

From our final data cube, we extracted line maps of Ha and 
[Nn]/16583, shown in Fig. [5^ and[5p. The Ha line map shows 
clumps which resemble the clumps of the F814W image. This 
is not surprising as both F814W(=rest frame UV) and Ha trace 
star formation. 



The [Nn]/I6583 line map shows a much more concentrated 
distributions, although [Nn]/i6583 is detected over the whole 
galaxy. If the inferred metallicity gradient is also present in the 
stellar population, this would put the highest metallicity and 
hence oldest stars in the centre of the galaxy, consistent with 
the redder color in the centre of the galaxy as found by Labbe et 
al. (120031 . The [Nn]^6583/Ha line ratio (0.2-0.4) is consistent 
with star formation over the whole galaxy. 



5. Velocity field 

We fitted a Gaussian to the Ha line at every position in the cube. 
We did this by hand using IRAF's splot because the S /N ratio 
of individual spectra is low and nearby OH line residuals con- 
fuse automatic fitting routines. For the outer parts of the galaxy, 
we were able to get a few additional data points by adding the 
spectra of 4 pixels (which were used as one pixel in the model- 
ing). The resulting velocity field is shown in Fig. [6^. It shows a 
remarkably regular spider diagram of a rotating disk. 



5. 1 . Tilted ring model 

We continued by fitting a simple model to the data using GIPSY 
(van der Hulst et al. [19921 Vogelaar & Terlouw 120011 which 
contains a software package designed for fitting tilted rings to HI 
velocity fields. Free parameters for each ring are in principle the 
centre of the galaxy, its inclination, the position angle of the ma- 
jor axis, systemic velocity, expansion velocity and the rotation 
velocity. For simplicity, we fixed the centre to the centre of the 
FIRES data, and the inclination and position angle of the major 
axis to those of the FIRES K band image. We also fixed the sys- 
temic velocity to the velocity of the central pixel, and fixed the 
expansion velocity to zero. This leaves only the rotation velocity 
as a free parameter for each ring. The resulting model is shown 
in Fig.|6p and the RC in Fig. [7] We will discuss the properties of 
the RC and the parameters of the the fit in the next section. 



5.2. Residuals from fit 

In Fig. |6j; we show the residuals of our tilted ring model fit. The 
lack of systematic residuals is a strong indication for a rotating 
disk. The mean absolute residual is 27 km s~ while large resid- 
uals (> 52 km s -1 ) only occur at the noisy outer edges of the 
galaxy. The mean residual excluding these points is 21 km s . 

These residuals are slightly larger than those of Wright et al. 
(2007) who observed a z ~ 1.5 disk galaxy and found a mean 
residual of 13 km s . They note that this is only slightly larger 
than the deviations observed for local spirals (~ lOkms -1 ). 
They explain their larger deviations by beam smearing. Our see- 
ing limited observation suffer even more from beam smearing 
than their AO observations. Modeling our residuals while taking 
into account these effects would require assumptions on the in- 
trinsic distribution of the Ha flux and dispersion and the intrinsic 
velocity field. We did not attempt to make such a model as the 
many uncertainties would not add much to the current conclu- 
sion: the small random residuals provide evidence that F257 is a 
disk galaxy with a regular velocity field quite similar to that of 
local galaxies. 

These results are in contrast to many other kinematic mea- 
surements of high-z galaxies which show significant deviations 
from regular rotation, e.g. the z = 2.38 galaxy for which Genzel 
et al. ( 120061 ) made a rotating disk model shows offsets up to 
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Fig. 5. a. Ha line map. b. [N ii]/16583 line map. Image size is again 3" times 4" and the images have been rotated so that east is to 
the top and north is to the right. The white crosses mark again the (K band) center (here determined from the PSF star positions). 




Fig. 6. a. The observed velocity field (not corrected for inclination). Some points at the edges are the average of four pixels, b. Best 
fit model velocity field, c. Residuals from the fit. d. Velocity dispersion field. The extent is smaller than for the velocity field due 
to the lower S IN. Instrumental resolution is 32 km s . Note that the velocity scale is different for this figure. The spatial scale and 
orientation is the same as in previous figures and the cross point marks again the (dynamical and photometric) center. 
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Fig. 7. Rotation curve from the velocity field model shown in 
Fig. |6j5. To guide the eye, the data points (with errorbars) are 
connected by the (black) thick line and the dotted lines are 
+ 10kms~ 1 and -lOkms 1 . The horizontal errorbars indicate 
the width of the rings. The full vertical errorbars are the errors 
on the fitted velocity, the dashed vertical errorbars are the range 
of velocities in the ring. The full (red) line is the best fit RC as 
discussed in Section l6~2l (a = 1.0, b = 0.5 and i = 54°). The 
long dashed (green) line shows an almost identical fit (a = 1.1, 
b = 0.5 and i = 64°). 



~ 170 km s 1 (comparable to the maximum rotation curve ve- 
locity) and deviations from a simple disk geometry. 

5.3. Dispersion field 

A peak in the measured velocity dispersion at the centre of ro- 
tation would be consistent with a rotating disk model. Although 
the S IN of individual spectra is poor, we made an attempt to 
measure the velocity dispersion field. We measured the disper- 
sion on the raw ID spectra, corrected for instrumental reso- 
lution and then smoothed the dispersion map with a 2 pixel 
width Gaussian to increase S IN. The dispersion map is shown 
in Fig [6}l. The dispersion is larger in the centre (maximum is 
103 km s ) and smaller in the outer parts (~ 65 km s ). The 
shape of the dispersion peak appears elongated, but the S /N is 
insufficient to make any statement about the precise shape. 

Peaks in the centre are not commonly observed at this red- 
shift: the z = 2.4 galaxy of Genzel et al. (2006 ) does show a peak 
in the centre but the overall shape is distorted, Forster-Schreiber 
et al. d20061 > find a peak for 1 out of 6 galaxies. Even at z = 0.6, 
Flores et al. (2006) find only 35% rotating disks (in their defini- 
tion: the major axis of the velocity field coincides with the major 
axis of photometry plus a peak in centre of the dispersion field) 
in a representative sample of emission line galaxies. 



5.4. Clumps 

The Ha line map and F814W image showed a clumpy structure. 
This raises the question to what extent the observed velocity field 
is biased by the inhomogeneous surface brightness distribution 
of Ha, emphasizing selected velocity ranges. The flux in clumps 
is 10-30 % of the total flux (measured in F814W), so although 
prominent in the image, they will not dominate the velocity field 
observed. We checked this in the model velocity fields where the 
F814W image was used as model for the true Ha surface bright- 
ness distribution (see Sect. 16. 2\ : there is no visible distortion of 
the velocity field at the location of the clumps. In addition, the 
smooth appearance of the velocity field and the agreement of 
kinematic and photometric (K band) major axis argues that the 
velocity field is dominated by the stellar mass. We conclude that 
the clumpy structure of the Ha surface brightness distribution 
has no significant impact on the velocity field observed. 

5.5. Summary 

The velocity field of F257 reveals a remarkably regular spider 
diagram. A tilted ring model, leaving only the rotation velocity 
as a free parameter for each tilted ring, fits the observations very 
well. The residuals from the fit are random and much smaller 
than in many other high-z galaxies with similar kinematic data 
(e.g. Genzel et al. 2006). The dispersion field shows an extended 
peak in the centre. The clumpy distribution of star formation ac- 
tivity has no significant impact on the velocity field observed. In 
the next section, we will discuss the RC found including the ef- 
fect of PSF smearing on Vfku and Rfuu, which we have ignored 
so far. 



6. Rotation curve 

The RC from the tilted ring model for the velocity field as dis- 
cussed in the previous section was shown in Fig. [7] It rises until 
it reaches its maximum velocity of ~ 170 km s (not corrected 
for inclination) at ~ 1" and then stays flat. Note that this veloc- 
ity is significantly larger than the velocity found in the synthetic 
long slit spectrum which is probably due to the integration over 
the width of the synthetic slit. The synthetic long slit spectrum 
(see section|4]and Fig. |4j» suggested that the RC flattens at both 
the receding and approaching ends at almost the same velocity. 
We repeated the fit of the velocity field as done in Sect. 15.11 but 
now fitting the approaching and receding sides seperately. The 
resulting RC, shown in Fig. [8] confirms that the RC flattens at 
the receding and approaching half at the same radius and ve- 
locity. Altough the Ha- flux distribution is not symmetric (see 
linemap in Fig. 0, this appears to have no significant influence 
on the best fit RC. For simplicity and S/N, we use the fit of the 
entire velocity field in the remainder of this paper. 

The errorbars on the RC were determined as follows: GIPSY 
calculates for each ring the error on the velocity AV and the 
range of velocities cry. The first is ~ 3 km s , the latter on av- 
erage 26 km s . Both AV and cry are shown in Fig. Q Here, 
equal weights are assigned to all pixels. Within a single ring, all 
pixels have approximately the same Ha flux and hence the same 
errorbar. The S /N drops from the centre outwards while at the 
same time the number of pixels in the fit increases. Assuming 
the errorbar on the velocity of a single pixel scales with the in- 
verse of the square root of the Ha flux, the errorbars on the RC 
are approximately constant except for the inner few rings and for 
r > l'/4 where the errorbars increase. 
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Fig. 8. Rotation curve from the velocity field model, approaching 
and receding side fitted separately. To guide the eye, the data 
points are connected by the thick line. The errorbars have the 
same meaning as in Fig. [7] 



6.1. Parameters of the fit 

We now turn to the parameters we fixed when we made a model 
of the velocity field: the position angle of the major axis, in- 
clination, centre, systemic velocity and expansion velocity. This 
is done in order to investigate if their values can be reproduced 
from the velocity field alone, and if so, if they are consistent with 
the values we assumed in Section lBTTl We also want to know the 
effect of small variations of these parameters on the best fit RC 
to get an estimate of the errorbars of the RC. It is not possible to 
leave all these parameters as free parameters at the same time. 
The number of free parameters is then too large for the number 
of independent data points, and the resulting best fit parameters 
for each ring vary strongly from ring to ring in an unphysical 
way (e.g. the inclination varies between 10 and 80°). Therefore, 
we study only the case of two free parameters: the rotation ve- 
locity plus one other parameter. 

1. If the angle of the major axis is left as the second free pa- 
rameter, one finds a consistent result between the orientation 
of the rings and the K band isophotes. This indicates that the 
orientation of the stellar mass distribution as traced by the K 
band light is consistent with the gravitational field as traced 
by the velocity field. Varying the position angle caused neg- 
ligible changes in the RC of the fit. 

2. When the expansion velocity is the second free parameter, 
we find V exp ~ km s~ with small scatter in all rings. 

3. Leaving the x and/or y position of the centre as free parame- 
ters), we find small (1-2 pixel) offsets from the K band cen- 
tre, varying from ring to ring and with no significant effect 
on the RC found. When we vary the fixed centre around the 
K band centre, we also do not find systematic changes in the 
RC found. We therefore keep the centre at the K band centre. 

4. We were not able to determine the inclination from the ve- 
locity field. If the inclination is a free parameter, the best fit 



inclination for each ring does not follow an ordered pattern, 
indicating that the S /N ratio does not allow determination of 
the inclination from the velocity field alone. However, fixing 
the inclination to different values causes small but system- 
atic offsets in the flattening velocity: more edge-on galaxies 
have slightly smaller flattening velocities than more face-on 
galaxies. The difference between i = 64° and i = 44° is 
~ 15 km s . Without the K band constraint on the inclina- 
tion, this small systematic difference would result in a larger 
uncertainty in Vf[ at . 

5. Leaving the systemic velocity as a free parameter, we find a 
small systematic offset from the velocity of the central pixel 
which we explain by the small difference in Vfi at of the re- 
ceding and approaching sides of ~ 10 km s _1 . The offset is 
smaller than the differences in Vfutt we find fitting with dif- 
ferent parameters (see conclusion below). Therefore, we ig- 
nore it in the remaining analysis. 

6. GIPSY also offers different weighting possibilities to the pix- 
els within a ring but this also does not significantly change 
the fitted RC. 

We were able to reproduce all parameters except the incli- 
nation from the velocity field alone. They all agree within the 
errorbars with the a priori assumend values. The agreement be- 
tween the centre and position angle of F257 in the Ks band image 
and the velocity field, and the agreement between the systemic 
velocity and the velocity of the centre, both strongly support the 
rotating disk interpretation. 

Changing the parameters in the fit except the inclination re- 
sults in random changes of the RC of the simplest model. The 
error on Vji at sin(/) derived by studying the variations due to the 
free parameter choice of the tilted ring model is ~ 15 km s . 
Varying the inclinations causes systematic offsets, that will be- 
come important when converting the observed velocities to true 
inclination corrected velocities. Therefore, we will focus on the 
inclination hereafter and keep all other parameters fixed. 

6.2. PSF 

So far, we have ignored smearing of our velocity field by see- 
ing. Although F257 is a large galaxy (we measure Ha over ~ 3" 
diameter), the seeing of ff.'55 has smeared the velocity field ob- 
served significantly. We will now model the effect of the PSF on 
the RC observed. For simplicity, we will focus on the RC only 
and not on the full velocity field. Specifically, we would like to 
know what the true V/uu sin(/) is, whether we can extract more 
information about the inclination of the galaxy and the radius 
where the RC flattens, Rji at . 

We generated our model velocity fields by taking the ob- 
served RC and modify it according to 



V mode i(r) 



a * V obs (br) if br < r„, 
a * V{r max ) if br > r m 



(D 



where V b s (r) is the velocity field observed at distance r from 
the centre, r max is the largest radius for we could measure the 
rotation velocity, V mo d e i is the velocity of the model, constant a 
scales the maximum velocity of the model and constant b scales 
the radius where the RC flattens. We then made synthetic model 
velocity fields (VF) for different a and b and created synthetic 
observations using 



VF, 



modeLobs 



(VF mo del,a,b X F, rue ) g PS F 

F trm ®PSF 



(2) 
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Fig. 9. The observed RC (open squares connected by the (black) 
dashed line) and the PSF corrected RC (full squares connected 
by the full (red) line). 



where VFmodeijobs is the model velocity field as "observed", 
VF mo dei,a,b is the intrinsic model velocity field for constants a 
and b and F true is the true flux distribution of the emission line. 
F trm is of course also unknown. We used the F814W image, as- 
suming that Ha and F814W (=rest frame UV) have the same 
distribution (as is likely, see section @). We used the elliptical 
gaussian as reconstructed from the PSF star images for the PSF. 
We then fitted a RC to each model for different inclinations. 

We found fitting RCs for a = 1 — 1.1. Larger values of a 
overpredicted V/uu sin(i) or R p a , (depending on the inclination), 
smaller values could not reproduce the maximum velocity ob- 
served. Our best fit Vfuu sin(/) is then 179 ± 9 km s~ . 

The best fit RCs have b = 0.5 - 0.6, although acceptable 
fits can also be found for slightly larger values up to b — 0.8. If 
Rfiat,obs ~ 1", Rfiat,true becomes 0.55 + 0.15" or 4.6 ± 1.3 kpc. 
Beam smearing has thus a profound effect on the observed flat- 
tening radius of the RC. This is illustrated in Fig. [9] where 
we show the observed RC and the PSF corrected RC. Altough 
Rfiat,true is of the order of the seeing, the radius where the RC 
flattens is resolved: in the velocity field, the (true) flattening ra- 
dius of the approaching and receding half are 2R fiat,true apart. 

Including the effect of beam smearing in our analysis, we 
were not able to constrain the inclination better as is shown in 
the almost identical fits with different inclination in Fig. [7] 

To summarise, we find Vp at sin(;) = 179 ± 17 km s , Rfim = 
4.6 +1.3 kpc and no better constraints on inclination after cor- 
rection for beam smearing. Using our best estimate of the in- 
clination (from the K band photometry) i = 54 + 5 °, we find 
V/iat = 221 ±22kms _1 . 

6.3. RC shape 

The RC of F257 flattens at 4.6 +1.3 kpc. In the local universe, 
the R fi at depends on the bulge-to-disk ratio, larger bulge-to-disk 
ratio's correlating with steeper rising RCs (Noordermeer 2006). 



The RCs of local early-type spiral galaxies rise to their maxi- 
mum velocity in < 1 kpc (Noordermeer 2006), which is clearly 
not the case for F257 as expected from the pure exponential 
disk profile found by Labbe et al. ( 12003I I. Swaters, Madore & 
Trewhella (2000) compare the RCs of their low surface bright- 
ness late-type galaxies with three late-type high surface bright- 
ness galaxies from Begeman ( 1987). They all flatten at 2 -5 kpc. 
So, the shape of the RC of F257 appears consistent with late type 
local spirals. 

Slowly rising RC at high redshift have also been noted 
by Swinbank et al. (2006). They studied gravitationally lensed 
galaxies at z ~ 1 using integral field spectroscopy and found that 
the galaxies with stable disk kinematics have slowly rising rota- 
tion curves and suggest that these have lower bulge-to-disk mass 
ratio's than their local counterparts. 



7. Masses 

From our model for the velocity field we calculate the total dy- 
namical mass of the galaxy within r = 1'.'5 



M dyn = 



(3) 



and find Md yn = 1.4 + 0.3xlO n M Q (depending on geometric 
assumptions, this equation is found with different constants of 
order 1 ; when we compared our finding with those of others, we 
used their constant.) 

Wuyts et al. d20071 l made an extensive study of stellar masses 
for high redshift galaxies in HDF-S and found for F257 3.5^[J \ x 
10 10 M (using a Salpeter IMF). Using a 'diet' Salpeter IMF, a 
Salpeter IMF with 30% less mass to correct for the overpredic- 
tion of low mass stars (as for example Bell & de Jong 120011 1. 
would reduce this mass to 2.5 x 10 M . 

A crude estimate of the gas mass may be found using the 
global Schmidt law for star forming galaxies which relates SFR 
and gas surface density (Kennicutt 119981 see also Bouche et 
al. 120071 for evidence that this relation also holds at high red- 
shift). Using the SFR calculated from the Ha flux and the to- 
tal area over which we observe Ha, this gives a total gas mass 
of 5.3 x 1O 1O M , and after correcting for extinction (Ay = 1.0, 
Wuyts et al. 120071 1 8.7 x 1O 1O M . Although the interpretation of 
this number is uncertain due to the large scatter in the Kennicutt 
relation of ±0.3 dex, it is comparable to or larger than the stellar 
mass, which is consistent with picture that this galaxy is building 
up its stellar mass at a fast rate. 

From these masses, we get the following mass ratio's 
(Salpeter IMF, extinction corrected): 

- M ga JM, = 2.5 

- MJM dyn = 0.25 ; 

- M baryon /M dy „ = (M» + M gas )IM dyn = 0.9 

Local galaxies with similar gas-to-stellar mass ratio are rare, and 
such galaxies have stellar masses ~ 1 - 2 orders of magnitude 
smaller than the stellar mass of F257 (Kannappan 2004). Similar 
gas-to-dynamical mass ratio's are locally only found in extreme 
late type galaxies (Binney & Merrifield 1998). The gas-to-stellar 
mass ratio is consistent with the z ~ 2 rest frame UV selected 
galaxies of Erb et al. (2006). The stellar-to-dynamical mass ratio 
is at the low end of the range found by Erb. et al. They find galax- 
ies with lower stellar-to-dynamical mass ratio's to be younger, 
and F257 is indeed a relatively young galaxy as found by Wuyts 
et al. Although F257 is photometrically young, its dynamical 
mass is at the high end of the range found by Erb et al. 
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The sum of stellar and baryonic mass is 1.2 x 1O U M which 
is almost the same as the dynamical mass. Within the errorbars, 
there is no evidence for dark matter within r < 12 kpc in this 
galaxy. 

F257 appears to be kinematically like a local late type spiral 
galaxy in many respects. We have seen that its velocity field is 
remarkably regular, a tilted model leaves small residuals and that 
the dispersion increases towards the centre. The RC slowly rises 
and then remains fiat like those of local late type spirals. F257 
is a very gas rich spiral compared to local galaxies. We will now 
further investigate the nature of this disk galaxy in terms of dy- 
namical support and specific angular momentum. We will then 
investigate if it follows the local B, K, stellar and baryonic mass 
TFR. 



8. Rotational support and specific angular 
momentum 

Galaxies can be supported by rotation and/or dispersion. Disk 
galaxies are supported by rotation. V/cr is a measure of the dy- 
namical support of a galaxy. There is no general agreement on 
what to use for cr for disk galaxies. Puech et al. ( 120071 1 use as 
cr the velocity dispersion in the z direction which they calculate 
from the mean observed cr in the outer regions of their disks. 
We followed Puech et al. (2007) and found an average cr z of 
57 km s and hence V/cr z = 3.9 ± 0.4. This is comparable to 
the median V/o~ z found by Puech et al., V/cr z = 3.8 + 2, for their 
Z = 0.6 rotating disk galaxies. Their local comparison sample 
has median V/cr, = 6.1 + 1.1 while their z = 0.6 perturbed rotat- 
ing disks have V/cr, = 2.4+2.5. Although the gas disk of F257 is 
heated compared to local galaxies, it is dynamically mainly sup- 
ported by rotation, comparable to the z ~ 0.6 galaxies of Puech 
et al. 00071 . 

An important question in galaxy formation theory is when 
disk galaxies acquired their angular momentum. The specific an- 
gular momentum of a thin rotating disk is 

y = 2R d V max (4) 

which gives log(- -2 — — ) = 3.2+0.2 for F257. This is slightly 

kpc km s 

smaller than the average specific angular momentum of galaxies 
of the same mass at redshift zero, but consistent with the specific 
angular momentum - velocity relation for local galaxies (Puech 
et al. [20071 Navarro & Steinmetz [2000t . Salucci et al. d2"007i l 
also found their z ~ 1 galaxies to be consistent with the local 
specific angular momentum - velocity relation. For other mass 
distributions, the factor 2 in equation |4] is slightly different, but 
this does not change our conclusion. 

Another measure of angular momentum is the spin parameter 
A . We followed again Puech et al. (2007 ) using their expression 
for A in terms of disk scale length and fiat rotation velocity and 
found A = 0.019+0.006 which is comparable to the values found 
by Puech et al. (2007 ) for their z ~ 0.6 rotating disk sample and 
their local reference sample and also close to the estimate for the 
local value made by Tonini et al. (2006) (A ~ 0.03). 

We conclude that F257 is a dynamically fully formed disk 
galaxy: it is mainly rotationally supported and has the same spe- 
cific angular momentum as local late type spirals. It is important 
to note that this makes F257 excellently suitable for a compari- 
son with the local TFRs, which we will do now. 



9. Tully-Fisher relations 

Having established that F257 is dynamically very much like a 
local disk galaxy, we compare it to local TFRs to compare its 
baryonic components to dynamically similar local galaxies. 

9. 1 . The rest frame B and K band TFR 

In Fig. [10] we plot the local rest frame B band and rest frame K 
band TFR of Verheijen (2001 ) together with our data point for 
F257. F257 is marginally brighter than expected from the local 
B band TFR. The offset is 0.7 + 0.4 mag in B band, where the 
errorbar reflects the uncertainty due to the velocity uncertainty. 
Correcting for extinction using the Ay found by Wuyts et al. 
(2007), the offset increases to 2.1 mag. For comparison: the lcr 
dispersion in the B band TFR as found by Verheijen d2001| l is 
0.41 measured along the magnitude axis. 

F257 is consistent with the local K band TFR with 0.4+0.4 + 
0.5 mag brightening where the first errorbar reflects the errors 
bars in the IRAC fluxes from which the rest frame K band mag- 
nitude was calculated (F257 is blended with a nearby source), 
and the second reflects the uncertainty due to the velocity uncer- 
tainty. For comparison: the lcr dispersion in the K band TFR as 
found by Verheijen d200Tb is 0.32. 

Local TFRs differ slightly from each other due to for ex- 
ample different choice of samples, different velocity parameters 
(e.g. linewidth or velocity, and also maximum or flat RC veloc- 
ity) and RC fitting methods (e.g. Puech et al. 2008 Kannappan & 
Barton |2004| l. The differences are smallest at the high mass end 
of the TFR and do not change our conclusions (van Starkenburg 
et al. 120061 and references therein). We choose the Verheijen 
(2001) TFR as reference TFR because it is based on velocity 
fields, it is defined in optical and near-IR bands and uses the flat 
RC velocity, which (locally) minimizes the scatter in the TFR. 

9.2. The stellar and baryonic mass TFR 

The relation between the stellar mass and rotation velocity is 
known as the stellar mass TFR. When we compare F257 with the 
local stellar mass TFRs of Bell & de Jong (2001) and Pizagno 
et al. (120051 ) as shown in Fig.[TTk we find that F257 lies below 
these stellar mass TFRs (after reducing the stellar mass found 
by Wuyts et al. 120071 by 30% to account for the different IMFs 
used). The offset in log(M«) is 0.66 + 0.20 and 0.44 + 0.12 re- 
spectively where the errorbars reflect the uncertainty due to the 
uncertainty in the velocity (the uncertainty in the stellar mass 
was 0.05). The lcr dispersion in the Pizagno et al. stellar mass 
TFR is 0.158 so the offset is significant at the ~ 2cr level. 

Wuyts et al. (2007) also discuss the effect on the best fit 
stellar mass of fitting a 2-component model instead of a 1- 
component model to the SED. The best fit stellar mass for a 
maximally old population plus a 100 Myr burst for F257 is 
1.0 x 10 11 M . This brings F257 in agreement with the local 
stellar mass TFR, as indicated by the arrow in Fig.QT] It is inter- 
esting to note that this 2-component model mass makes F257 a 
nearly (~ 75%) maximum disk. 

Some authors (e.g. Verheijen I200T1 Bell & de Jong l200Tb 
added the gas mass to the stellar mass to get a baryonic mass 
TFR, see also Fig. [TTb . Although the gas mass of F257 is un- 
certain due to the large scatter in the Kennicutt relation, we also 
compare F257 with these local baryonic mass TFRs and find that 
F257 is consistent with the baryonic mass TFR from Bell & de 
Jong (1200Tb : log(M /w ) = -0.04+0.29+0.3 where the first error- 
bar reflects the uncertainty due to the uncertainty in the velocity 
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Fig. 10. The rest frame B band (left) and rest frame K band (right) TFR. The arrow indicates the extinction correction using the best 
fit Ay from Wuyts et al. (2007). The local TFR and its scatter (full and dashed lines) are from Verheijen (1200 11 1. The errorbar on the 
absolute B magnitude is smaller than the point size and is therefore not shown. 



and the second errorbar reflects the uncertainty in the baryonic 
mass (assuming of a factor 2 uncertainty in the gas mass). 

Although Bell & de Jong ( 1200 11 1 use an earlier version of 
the data set of Verheijen (2001), they find significantly different 
baryonic mass TFRs, especially at the high mass end. This differ- 
ence is explained by the different M/L that they use to convert K 
band light to stellar mass. The (constant) M/L used by Verheijen 
is ~ 3 times larger than the (best fit for each galaxy) M/L of 
Bell & de Jong (1200 11 1. Verheijen (1200 11 1 uses this M/L because 
it minimizes the scatter in the baryonic mass TFR. However, ac- 
cording to Bell & de Jong, this M/L is unrealistically high. A 
lower M/L would shift the Verheijen relation downward, espe- 
cially at the high mass (i.e. low gas fraction) end. 

9.3. Comparison to other high-z samples 

Several authors have attempted to measure rest frame near-IR 
and stellar mass TFRs at high redshift. Conselice et al. ([2005 ) 
analyzed 101 disk galaxies in the redshift range 0.2 < z < 1.2 
and find no evolution in the rest frame K band TFR and in the 
stellar mass TFR. 

Kassin et al. ( 120071) use a combination of rotation velocities 
and velocity dispersions (S 0.5 = ^jo.5V^: ot + cr^ as ) to construct a 

stellar mass TFR for galaxies with strong emission lines in the 
redshift range z = 0.1 — 1.2. This reduces the scatter in the stel- 
lar mass TFR. However, they find no z-evolution in the stellar 
mass TFR using this S 0.5 parameter. F257 is a rotation domi- 
nated galaxy and hence consistent with their normal stellar mass 
TFR and their S 0.5 parameter stellar mass TFR. 

Weiner et al. (|2006) study the rest frame B and J band TFR 
out to z — 1.2. They find significant brightening in rest frame 
B band (1-1.5 mag between redshift 0.4 and 1.2). However, al- 
though they find that the slope of the J band TFR evolves, the 
overall luminosity in J band does not. 



Flores et al. (|2006) published the first high redshift TFR us- 
ing integral field spectroscopy. They observed a sample of [0 11] 
emitting galaxies at z = 0.4 - 0.75 using the GIRAFFE IFU on 
the VLT. They find that only 35% of their sample are rotating 
disks in their definition (the major axis of the velocity field coin- 
cides with the major axis of photometry plus a peak in centre of 
the dispersion field). They find no evolution and small scatter in 
the rest frame K band and stellar mass TFR and some galaxies 
that are brighter in rest frame B band. 

To summarize, other authors have not found evolution in the 
near-IR or the stellar mass TFR out to redshift 1.2. This non- 
evolution in near-IR z — 1 TFR was predicted by Buchalter, 
Jimenez & Kamionkowski (2001) and in the stellar mass z = 1 
TFR by Portinari & Sommer-Larsen (2007). Our findings sug- 
gest that the K band TFR may not evolve out to redshift 2. F257 
is ~ 2cr below the stellar mass TFR, a hint that the stellar mass 
of F257 is still building up, which would be consistent with its 
young age, large SFR and gas mass. However, it should be noted 
that the offset depends on the SED modeling. If a maximally 
old population is added, F257 is consistent with the local stellar 
mass TFR. 



9.4. Summary and discussion 

We conclude that F257 is a fully formed, rotationally supported 
disk galaxy. Within the errorbars, it is consistent with the local 
K band TFR. The small offset from the local B band TFR can be 
explained by star formation activity and the young stellar popu- 
lation. 

F257 lies slightly below the local stellar mass TFR, an off- 
set significant at the 2<x level. Evolving from redshift 2 to the 
present, the stellar mass must increase without a significant in- 
crease in the K band luminosity. However, the best fit stellar 
mass strongly depends on the SED modeling. When we add a 
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Fig. 11. a. The local stellar mass TFR and lcr scatter of Pizagno et al. (2005} (long dashed and dotted lines) and Bell & de Jong 
(2001 ) (dashed line), b. The local baryonic mass TFR and lcr scatter of Verheijen (2001) and Bell & de Jong (2001 ) (dashed line). 
The arrows indicate where F257 would lie for the maximally old plus young burst assumption (see text). 



maximally old population, F257 is in agreement with the local 
stellar mass TFR. 

F257 is consistent with the baryonic mass TFR, but the un- 
certainties are large. To study offsets from the baryonic TFR, 
one needs more galaxies and/or more accurate gas masses. We 
present our full sample in van Starkenburg et al. (2008b). In the 
near future, ALMA will provide us with the latter. 

We finally note that the regular dynamcial sctructure of F257 
is remarkable for its young age: the estimated age of the stellar 
population of ~ \6QMyr equals the rotation period of the galaxy 
at 1 - 2Rd and at 1 - 2Rfi al lrue ). Therefore, the galaxy has dy- 
namically relaxed on a very short timescale, much shorter than 
the timescale for the buildup of the final stellar population. 

10. Summary and conclusions 

Labbe et al. (2003 ) presented photometric evidence that F257 is 
a large disk galaxy at z = 2.03. We confirm the disk nature of 
this galaxy by its kinematical properties and find that is in many 
respects similar to local disk galaxies. Specifically, we find: 

1 . The velocity field is consistent with a rotating disk galaxy, a 
tilted ring model leaves small random residuals and no evi- 
dence for non circular motions; 

2. There is no evidence for an AGN in the form of broad lines 
or a large [N n]/Ha ratio; 

3. The [N n]/Ha ratio peaks in the centre; 

4. The disk is primarily rotation supported: V/cr, ~ 3.9; 

5. The dispersion increases towards the center; 

6. The RC slowly rises like that of a late-type spiral in the lo- 
cal universe and then flattens. There is no evidence for dark 
matter within r < 12kpc; 

7. F257 has a specific angular momentum and spin parameter 
comparable to those of local galaxies; 

8. F257 is consistent with the local K band TFR; 



9. F257 lies below the local stellar mass TFR at the 2<x level 
and is consistent with the local baryonic mass TFR. A 2- 
component model for the SFH increases the stellar mass so 
that F257 is consistent with the stellar and baryonic mass 
TFR and makes F257 a nearly maximum disk. 

Despite its very young age, high star formation rate, gas mass 
and clumpy star formation, F257 is very similar to local late type 
spiral galaxies. Its dynamical properties are more like those of 
local galaxies than those of other galaxies observed at similar 
redshift. 

These findings put constraints on the evolution of F257 from 
Z = 2.03 to the present. If the stellar mass increases so that F257 
becomes consistent with the local stellar mass TFR, the K band 
mass-to-light ratio must also increase so that that the total K lu- 
minosity does not change. Infalling gas must have the same spe- 
cific angular momentum as the galaxy, i.e. it must fall in from 
the outer parts of the present halo. 

Acknowledgements. LvS wishes to thank Juha Reunanen for lengthy discussions 
about SINFONI data reduction. We thank Eric Bell for an interesting discussion 
and Marc Verheijen for his help with the local TFRs. 
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Appendix A: FOV reconstruction 

A. 1. Introduction 

Our targets are too faint to acquire directly or to detect in individ- 
ual frames. When we made our first observations with SINFONI, 



32 




X 

Fig. A.l. SINFONI observing strategy. The target is observed at 
positions 1, 2, 3 and 4. The distance between position 1 and 3 is 
I", between 1 and 4 4". The positions of the slitlets and the X 
and Y directions are marked. Slitlets are numbered from top to 
bottom. The total FOV size is 8"times 8". 



we adopted the normal strategy for these faint targets: blind off- 
setting from a nearby bright star. We also followed the recom- 
mendation in the SINFONI manual to reobserve the blind offset 
star in the same OB as the science observations to be able to ac- 
curately combine the data from different OBs. We refer to these 
star observations as "PSF stars" (the template used for the obser- 
vation is the calibration template "PSF"). 

We also followed the standard procedure to reconstruct the 
FOV from the raw images. SINFONI's FOV is sliced into 32 
slitlets (see also Fig. IA.lt that are combined into a single pseudo- 
slit and then dispersed. These 32 slitlets are projected on the 
detector in a brick-wall pattern. The slitlet distances (measured 
on calibration data made with fibers) define the distances on the 
detector between points with the same x coordinate (see Fig. IA.ll 
for a definition of the x and y direction). These slitlet distances 
are used to reconstruct the FOV in the reduction process. 

Having reduced our first data sets, we noticed that the re- 
constructed images of the PSF and standard stars showed some 
peculiarities. The centre of the star changed from slitlet to slitlet, 
often back and forth, but in some cases the gradual shift in cen- 
tre resulted in a stretched stellar image. An example is shown in 
Fig. IA.2I Also, the overall shape of the stars was not round but 
elliptical, especially if the seeing was good. FWHM ratios var- 
ied between 1 and 1 .6, much larger than we could explain. As all 
reconstructed stars looked different, a simple error in the slitlet 
distances was excluded. 

After careful study of the reconstructed images of our PSF 
and standard stars of our first observing runs, we concluded that 
the details of SINFONI's FOV reconstruction are more compli- 
cated and that the quality of the velocity fields measured might 
be improved if the details are treated carefully. In this Appendix, 
we describe our investigations, the effects we found, their effect 
on the data quality and our new improved observing strategy. 
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Fig. A.2. The reconstructed image of one of the bright standard 
stars observed for F257. To guide the eye, the center of the star 
in each slitlet has been marked. 



A.2. Further investigations: varying pixel scale 

To investigate the details of SINFONI's FOV reconstruction, a 
PSF star was observed at the same positions in the FOV as the 
science target, in some cases even before and after the science 
target. The positions in the SINFONI FOV are shown in Fig. 
|A~T1 

When we compared the star traces of the PSF stars that are 
in the same slitlets in the raw frames (i.e. stars 1 and 3, and 
stars 2 and 4 in Fig. IA.lt , we noted that the distance between the 
star traces varies as is illustrated in Fig. IA.3I From the telescope 
offset in the x direction (1") and SINFONI's spatial pixel scale 
((X'125 pix~') the distance between the star traces should be 8 
pixels in all slitlets. The average value we found was indeed 8 
pixels, but we found offsets of ~ 10% around this value (see Fig. 
|A.3t . In other words: the pixel scale varied up to 10%. The effect 
was not stable in time, we have seen significant changes in only 
1 hour as can also been seen in Fig. IA.3I 

We emphasize that we measured these positions on the raw 
frames, so that reduction artifacts are excluded. Also, this effect 
is not caused by errors in the telescope offset: the pixel scale also 
varied (by less than 10%) from slitlet to slitlet for a single pair 
(see the upper or lower half of the panels in Fig. |A.3I >, and the 
pixel scale in the y direction (i.e. the distance between stars 1 
and 4 or 2 and 3 in Fig. lA.U was stable to measurable accuracy. 

This explains the odd shapes of the reconstructed stellar im- 
ages. The slitlet distances measured on fiber data are only valid 
for that observation and only approximately valid for other ob- 
servations. The changing pixel scale explains why the FWHM 
of the stars in x and y direction is not the same for good seeing 
data. 

SINFONI is a Cassegrain-mounted instrument, which makes 
is sensitive to flexure. As we saw in later observing runs (where 
we adopted our new multiple PSF star strategy (see below)) in 
many consecutive hours of observations of the same target, the 
pixel scale changes smoothly in time and hence parallactic angle. 
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Fig. A.3. The star traces for the test observing sequence with ob- 
serving strategy as shown in Fig. IA.ll Each panel represents one 
slitlet, the slitlet number is indicated in the lower left of each 
panel. The vertical lines show the positions of the star traces be- 
fore (full lines) and after (dashed) lines observations of the sci- 
ence frames (time difference is about 1 hour). For each slitlet, the 
position of the first observation of the left star has been shifted to 
zero. The centers of the two stars could be measured in 8 slitlets. 
Those of the two stars in the lower half of the FOV (stars 1 and 
3 in Fig. lA.U are in the upper 8 panels, those of the stars in the 
upper half of the FOV (stars 2 and 4 in Fig. lA.U are in the lower 
8 panels. The distance between two star traces should be 8 pixels 
if the pixel scale is (X'125. Note that the distances between the 
two stars are different and change differently in each slitlet. Also 
note that the distance is less than 8 pixels for most slitlets in the 
lower half of the FOV and more than 8 pixels for most slitlets in 
the upper half of the FOV. 



Therefore, we suspect flexure in the instrument as the origin of 
the effects described above. This could also explain the 'moving 
slitlets' phenomenon: for some parallactic angles, the positions 
of the slitlets on the detector shifts up to ~ 2 pixels in a few 
hours (of observing the same target). 

We do not have sufficient observations of PSF stars to be 
able to predict the pixel scale for a given slitlet at a given tar- 
get at a given time (or even to be sure that pixel scale can be 
predicted from parameters like parallactic angle and airmass). 
We have indications that the upper and lower half of the FOV 
have systematically different pixel scales (as is the case for the 
example shown in Fig. IA.31 >. We only looked at data taken with 
SINFONI's largest FOV without AO and have not searched for 
differences between J, H, H+K and K band although we have 
noted these effects in all 4 filters. 



A.3. Effects on data quality 

Given these complex pixel scale changes, the question arose 
what the consequences were for data sets with observation of 
a single PSF star and what improvements could be made in fu- 
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ture observations (including the observations discussed in this 
paper). 

We assumed that the addition of many frames will average 
out the pixel scale differences so that the pixel scale in the com- 
bined image is the same in both spatial directions: 07125 (which 
seems safe to assume, as the pixel scales we find have on average 
this value). The changes in object size are negligible compared 
to the seeing even for the large high redshift galaxies in our sam- 
ple, moreover seeing variations are usually larger than 10%. 

However, the varying pixel scale results in a position uncer- 
tainty which can be important. In our science exposures, the tele- 
scope offset between A and B frames in the x direction is about 
4" for all our data (this paper, and van Starkenburg et al. 2008a 
2008b). 10% of 4" is comparable to our best seeing observations 
(e.g. this paper 0755). This position uncertainty causes an addi- 
tional smear in the x direction which cannot be removed with 
the observation of one PSF star per OB. As our targets are too 
faint to measure their position in individual frames, our new ob- 
serving strategy should solve this position uncertainty in another 
way. 

To summarize: the varying pixel scale causes a position un- 
certainty that could degrade the data quality in one spatial direc- 
tion if the target is faint and observed under good seeing condi- 
tions. 



tortion entirely. As we are only interested in a single emission 
line, we solved this by measuring the position and shape of the 
PSF stars at the wavelength of the emission line. 



A.4. Multiple PSF star observing strategy 

Given these results, we changed our observing strategy 
as follows. The position uncertainty is the most impor- 
tant factor to correct. Observing a PSF star at every 
position where we observed the science target and as 
closely in time as possible provides us with accurate po- 
sitions (i.e. in term of an ABB'A'-like observing cycle 
ApsfBpsf ^sciencelJscience A sc j ence B sc j ence A p S fB p S f). 

We then use these PSF positions when we combine our data. 

As a bonus, we can seeing-weight the images (using the 
FWHM in the y direction) to improve image quality further. 
Another benefit is that we do not have to worry about the 'mov- 
ing slitlet' effect anymore (it is unclear whether the moving 
slitlets are simply a shift in projection of the same part of the 
sky or a projection of a slightly different part of the sky or a 
combination of both). 

When preparing the observations, we also take into account 
that the FOV size might be slightly smaller in the x direction than 
expected or that the target might be closer to the edge of the FOV 
than you would expect from the offsets. Also, when possible we 
put the major axis of the galaxy along the y direction in the FOV 
to ensure maxium resolution where we need it most. 

We did not attempt to correct the individual frames and 
slitlets for the pixel scale variations as we saw before that these 
effects are negligible compared to the seeing. Moreover, in the 
case of correcting individual slitlets, one would need very bright 
PSF stars to get enough flux in all slitlets where the (large) sci- 
ence target is, and those are in general not available (the example 
shown in Fig. |A.3| is with a very bright PSF star). 



A.5. Other artifacts 

For the very bright standard stars, we noted that their centres 
in the reconstructed image varied as a function of wavelength. 
This may indicate even more complex flexure effects, but a more 
probable explanation is that the distortion correction of the (ear- 
lier versions of the) SINFONI pipeline did not remove the dis- 



